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ABSTRACT 

The strong degeneracy of the ^^C ignition layer on an accreting neutron star results in a hydro- 
dynamic thermonuclear runaway, in which the nuclear heating time becomes shorter than the local 
dynamical time. We model the resulting combustion wave during these superbursts as an upward 
propagating detonation. We solve the reactive fluid flow and show that the detonation propagates 
through the deepest layers of fuel and drives a shock wave that steepens as it travels upward into lower 
density material. The shock is sufficiently strong upon reaching the freshly accreted H/He layer that 
it triggers unstable "^He burning if the superburst occurs during the latter half of the regular Type I 
bursting cycle; this is likely the origin of the bright Type I precursor bursts observed at the onset of 
superbursts. The cooling of the outermost shock-heated layers produces a bright, w 0.1 s, flash that 
precedes the Type I burst by a few seconds; this may be the origin of the spike seen at the burst onset 
in 4U 1820-30 and 4U 1636-54, the only two bursts observed with RXTE at high time resolution. The 
dominant products of the ^^G detonation are ^*Si, '^^S, and ^^^Ar. Gupta et al. showed that a crust 
composed of such intermediate mass elements has a larger heat flux than one composed of iron-peak 
elements and helps bring the superburst ignition depth into better agreement with values inferred 
from observations. 

Subject headings: accretion, accretion disks — nuclear reactions, nucleosynthesis, abundances — stars: 
neutron — X-rays: bursts 



1. INTRODUCTION 

Superbursts are powered by unstable thermonuclear 
burning on the surface of an accreting neutron star in 
a low mass X-ray binary (for reviews, se e lKuulkersll200l 



ICummingj[200HStrohmaver fc Bildstenll2006f) . Their re- 
currence time (~ 1 yr) and duration (pa hours) suggest 
ign ition densities of 10* — 10^ g cm~^ (jCumming et al.l 
[200^ 1. 100 - 1000 times larger than that of nor- 
mal Type I X-ray bursts. If t he superburst fuel is 
^!C (|Cumming fc BildstenI l2001t iStrohmaver fc Brownl 
l2002t l. then the observed energy release (w 10^^ ergs) 
implies a ^^C mass fraction X1 2 > 0.1 in the m atter ac- 
cumulated from Type I bursts (jCumming et a l. 2006). 

Previous studies assumed the entire layer of ^^C 
fuel burns instantly and hydrostatically, and ob- 
tained a good match to the late-time (> 10'^ s) 
light curves and 
normal bursts 



accounted for the quenching of 
t he days following a superburst 



(ICumming fc Bildsteiill2001l: [Strohmaver fc Brownl 



Gumming fc MacbethI 12001 ICooper fc Naravanf 



2002 



2005 



Gumming et al .1 120061 ). They did not, however, flt the 



early-time light curves nor provide a mechanism to trig- 
ger the Type I bursts that precede superbursts by « 
10 s. Such precursor bursts have been found in all 
five cases in which t he onset of the superburst was 
observed (4U1820-30: IStrohmaver fc Brownl 120021 : 4U 
1636-54: [Stroh maver fc Markw ardtl 120021: KS 1731-26 : 
iKuulkers et all 2002: 4U 1254-69: lin't Zand et al.ll2003l : 
GX 17-1-2: li~Zand et al.ll2004bl) . 



In Weinberg, Bildsten, fc Brown (2006; hereafter 
WBB), we showed that the strong degeneracy of the su- 
perburst ignition layer results in a hydrodynamic ther- 
monuclear runaway, in which the heating time th = 
(dlnTb/dt)"^ becomes shorter than the local dynami- 
cal time td — h/cs ~ 10~^ s, where Tb, h, and Cg are 
the temperature, pressure scale height, and sound speed 
at the base of the burning layer. We start in § |2| by 
describing the conditions for the spontaneous initiation 
of a detonation, and show that the observed superbursts 
are likely deep enough into the neutron star that their 
plane-parallel burning will initiate a detonation. We thus 
model the hydrodynamic combustion wave that forms 
once th < as a detonation, showing its impact on the 
overlying surface layers. 

We show (§ |2|) that as the detonation propagates up- 
ward, it drives an outgoing shock wave that steepens as 
it travels into lower density material. When the shock 
impacts the freshly accreted H/He layer, its overpressure 
is Ap/p K, 10. Our calculations (§|2|) suggest that about 
half the time there is enough H/He present in the layer 
that the shock can ignite the ^He and trigger a Type I 
precursor burst. We compute full light curves in §|4|and 
show that a detonation with a shock-triggered Type I 
burst can explain features of the early-time light curves 
while still accounting for the late-time behavior. We also 
show that the cooling of the outermost shock-heated lay- 
ers results in a bright, sub-second, flash that precedes 
the precursor by a few seconds. This could be the origin 
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of the spike seen at the burst onset in 4U1820-30 and 
4U 1636-54, the only two bursts observed with RXTE at 
high time resolution. 

We show in § |4] that even if the shock does not ignite 
the ''He, the cooling of the ashes of ^^C burning at depths 
just below the H/He layer can result in a precursor that 
appears similar to a Type I burst, but with a peak lumi- 
nosity that is smaller by a factor of ~ 2 — 3. This may 
explain why three out of the five precursor bursts were 
only aboTit half as brig ht as the ordinary Type I bursts 
pCuulk ers 2004 : Gumming et al.1 .2006) . 

2. HYDRODYNAMIG CARBON BURNING 

After the thermally unstable ignition of ^^G, there is an 
hour-long convective stage where the energy generation 
rate from ^^G burning, ec, is sufficiently small that the 
heating time is much longer than the eddy turnover 
time te = h/vc, where Vc = Vc(Tb) is the typical velocity 
of a convective cell (WBB). For low ^^G ignition depths 
(i.e., low densities) or initial ^^G mass fractions X12, the 
^^G completely burns in this convective stage, and hy- 
drostatics describes all stages of burnii ig (just as with 
'^He burning during Type I X-ray bursts; IWeinberg et al.l 
HoOGb). 

However, for larger ^^G ignition depths 2/6 ^ 2 x 
10^^ g cm~^ (ignition densities^ pb > 1.5 x 10* g cm^'^) 
and X12 > 0.2, ec becomes so large that th < te and 
convection becomes inefficient at transporting the en- 
ergy release outward (uc — O.OTc^ when th — t^ for 
yf, = 10^2 g cm-2, X12 = 0.2; WBB). As a result, the 
fuel is rapidly consumed within a thin layer near the base 
in a local thermonuclear runaway, during which time th 
becomes shorter than the dynamical time td and hydro- 
dynamic burning begins. The transition from th = te to 
th = td only requires a w 20% increase in temperature 
(WBB). We calculate the minimum heating time th,min 
using the estimate in WBB (equation [5]) and in Figure 
[T]plot the relation between yj, and X12 where ih,min = td- 
Ignitions above this line will become hydrodynamic. 

The only direct information we have about yi, and X12 
is that derived by Gumming et al. (2006) from fits to six 
observed superburst light curves. These are plotted in 
Figure [T] assuming a 30% uncertainty in the distance to 
each source (corresponding to the uncertainty in the dis- 
tance to 4U 1254-690, the only one of the six sources 
with a reported best fit distance with error bars).^ The 
measured values of yi, and X12 suggest that the ^^G likely 
burns hydrodynamically. We argue in § 12.11 that a possi- 
ble (but by no means proven) outcome is a spontaneous 
initiation that leads to a propagating detonation. 

2.1. Spontaneous Initiation of a Detonation 

Even for chemical combustion in gases, it re- 
mains challenging to theoretically predict the condi- 
tions needed for a spontaneous detonation (jLed I1977L 

^ Assuming a relativistic degenerate electron gas, the density at 
the base is pi, ^ (5 X 10* g cm-3){gi4/2)3/4{0.5/ye)s/?'Cl where 



Ignition Density (10^ g cm ^) 



514 = a/iQi 



is the gravity, Yc, is the electron fraction, and 



2/6,12 = Vb/W^^ g cm 2, 

The light curves constrain the nuclear energy release £17 = 
ergs g~^. To estimate X12 we assume X12 = O.I-E17 
corresponding to ^^C burning to ^^Fe. As we show in ^ 12.31 how- 
ever, the ^^C may only burn to ^*Si, in which case X12 is larger 
by f5! 40%. 
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Fig. 1. — Range of superburst ignition parameters yi, and X12, 
with the top axis showing the base density at ignition, pi,. The 
region above the solid line has a minimum heating time th.min l^ss 
th an the dynam i cal tim e t^. The six data points are from fits by 
G umming et al.1 I I2006I ) to observed superburst light curves. The 
error bars assume a 30% uncertainty in distan ce to each sourc e 
(and use the scalings with distance given in Gu mming et al.ll2006i ') . 
The open star at yi, = 10^^ g cm~^, X12 = 0.2, indicates the 
parameters chosen for the detonation calculation shown in Figures 
2-5. The upper (lower) dashed curve shows where Kmax'i = 5 

(^max^ — !)■ 



11984 iBdzil fc Stewart! 120071 ). However, a widely ac- 
cepted criteria for initially addressing the question 
is the requirement that a sufficiently large volume 
(defined below) burn in le ss than a sound crossing 
time (B linnikov & Khokhl ovl Il986l. [l987t IWooslevI 119901: 
iNiemever fc Wooslevlll997i : iKhokhlov et allligm 

We first estimate the size of the region that burns 
in less than a sound c rossing time. According to the 
Zel'dovich criterion (see lBartenev fc Gelfandl[2000l for an 
updated discussion), such burning requires a spatial re- 
gion where the temperature gradient is so shallow that 
nearly simultaneous runaways can occur. The quanti- 
tative criterion is that the phase speed of the burning 
Vph = \dtnuc/dx\^'^ > Cs, where t^^c = (dlnec /dt)~^ = 
th/v ~ CpT/vtc is the time scale for the accelera- 
tion of nuclear burning and v = dlnec/dlnT « 21 at 
T = 2 X lO^K. We first evaluate this in the vertical direc- 
tion assuming an adiabatic temperature gradient at the 
moment when tnuc = te, which occurs at ?» 2 x 10^ K 
(we integrate over the entire convective region when com- 
puting inuc; see Figure 1 and equation [3] in WBB). 
Since the heat capacity is determined mainly by the elec- 
trons, Cp oc (|Gumming fc BildstenI [2OOII ) . and 
thus inuc Of T'^^^ jy. Requiring that Wph > Cg then im- 
phes 

th{y) < Xtd « 3td, (1) 

where A = + {1^ — 2)Vad] and Vad = 

((ilnr/(iln?/)ad ~ 0.3 is the adiabatic index. At run- 
away th{yb) — td, and Uph > Cs between the base and a 
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height 



Zcrit ~ /iln(j/f,/ycrit) 



AlnA 



0.2/1. 



(2) 



Above that location, the temperature is not adequate to 
help initiate the detonation. Of course, because of the 
small scale height compared to the radius, there is a much 
smaller temperature gradient in the horizontal direction, 
an issue we will address momentarily. 

If the initiating region calculated above is too small, 
the initial shock is unable to trigger further burning be- 
fore it decays by geometrical dilution. Direct numerical 
simulations and experiments always find that to trigger 
a detonation the initiating region must be a few orders of 
magnitude longer than th e induction length of the deto- 
nation lind ~ vcjtind (see lEck ett et al.l!200Cl| for a recent 
discussion). Here vcj is th e Chapman- Jouguet veloc- 
ity (see e.g., lKhokhlovlll989l) and tjnd is the induction 
time, the time scale to complete the burning behind the 
shoc k. In the case of a spherically symmetric ignition 
(see iHe fc Clav in 1994), such as in a centrally ignited 
Type la supernova, this minimum distance is ex pressed 
in terms of a maximum radiu s o f curvature Atmax j Sharpd 
[2OOI iDursi fc Timmesl[2006l) . iDursi fc TimmesI (|2OO60 
carried out direct numerical simulations for carbon det- 
onations, and derived a fitting function for Kmax(p, -^^12)- 
The most conservative estimate of the critical conditions 
needed for a detonation is to then demand that the ver- 
tical length of the spontaneously igniting region, Zcrit, be 
larger than k^\^. By equation ([2]), this translates into a 
requirement that K^ax^ > 5, as plotted in Figure [T] Most 
of the data lie below this line, making it challenging to 
trigger a spherical detonation by a vertical initiation. 

Building on the arguments for Type la supernovae core 
detonations, our discussion has assumed that the vertical 
temperature gradient (where the temperature decreases 
by « 30 — 40% over a pressure scale height h) plays the 
critical role. However, the superburst occurs in a thin 
shell at the surface of the neutron star, and the hour 
long convective phase prior to the hydrodynamic run- 
away should be adequate to establish transverse temper- 
ature gradients smaller than the vertical gradient. Only 
a numerical anelastic calculation could accurately calcu- 
late such a gradient, but we think it reasonable to assume 
that the initiation length is larger in the plane, poten- 
tially exceeding h at some location. Such a cylindrical ge- 
ometry also has the advantage of less severe geometrical 
dilution, where all analyti c calculations yield a factor o f 
two smaller value for «;~J^^ (|Leelll977HHe fc Clavinlll994D . 
For this reason, we plot an additional line at Kmax^ = 1 
in Figure [TJ representing our current best guess of the 
constraint for detonation initiation in a planar geometry. 
All the data points are above this line, as is the explosion 
we are simulating in this paper (shown by the open star). 
Although more work clearly remains to fully assess the 
onset of a detonation in a planar geometry after a long 
convective phase, we will proceed assuming that it can 
happen. 

2.2. Numerical Method 

We assume a plane parallel geometry with constant 
gravity g — 2.4 x 10^^ cm s^^ (corresponding to M = 
I.4M0, i? = 10 km) and treat the detonation in only 
the vertical direction, parametrized by the column depth 



y into the star (units of g cm~^), where dy = —pdr. 
Modeling the lateral propagation is imp ortant and left 
for future study fsee lZingale et al.ll2001l who consider a 
laterally propagating ^He detonation in the ocean of a 
neutron star). We solve the reactive fluid flow equations 
using a one-dimensional, explicit, La grangian, finite dif- 
ference scheme (see e.g., iBen j IT99ll) . We use operator 
splitting to couple the hydrodynamics to a nuclear en- 
ergy generation network. The network contains 13 iso- 
topes: {4He, ^^C, i^o, 20Ne, ^^Mg, ^^Si, 32s, 36 Ar, ^oca, 
44rp|^ "^^Cr, ^^Fe, ^^Ni}, and includes a-chain, heavy-ion, 
and (a,p)(p, 7) reactions^. We assume electrons, ions, 
and photons supply the pressure and calculate the equa- 
tion of state, volu metric neutrino emissivity, and thermal 
conductivity, as in iBrownl (|2004f ). 

We model the detonation and shocks over the col- 
umn depths 10^ — 10"'^^ g cm~2 using 4000 grid points 
spaced uniformly in logy. The ^'^G layer spans the 
range j/hc < y < yt ^ 10^^ g cm"^. The depth of 
the base of the freshly accreted H/He layer, ytici'P) = 
0yHc(0 = 1) ~ 10^ g cm~2^ varies linearly with the phase 
(f) — i/irecur of the Type I burst cycle, whose recurrence 
time irocur 10'' — 10^ s. We assume the accreted gas has 
a solar composition and consider local accretion rates per 
unit area in the range rh = 0.05 — 0.3TOEdd, where riiEdd 
is the local Eddington accretion rate (i.e., rough ly the 
range inferred from observations: iKuulkerd 120041). The 
^^C layer has a mass fraction 1 — X12 of ^^Fe and for 
y > yb the composition is pure ^^Fe. For a given y^, X12, 
(j), and TO, we solve for the crustal heat flux that results 
in unstable ignition at yb (see WBB) and then c alculate 
J/Ho(0 = 1, 'Ti) as in lCumming fc BildsteiJ (j2000D . 

The shock reaches the top of our grid {y = lO'* g cm^^) 
at a time itop ~ 10 /iS after the onset of the hydrody- 
namic runaway. The underlying layers are adiabatically 
expanding at t = /top, though they remain bound to the 
neutron star and in plane parallel (see § 12. 4p . We assume 
that for t > ttop the flow remains adiabatic and that after 
a few dynamical times (tens of ^s), the shocked layers set- 
tle into a new, puffed-up, hydrostatic equilibrium, with 
an entropy profile given by that at t = itop- This post- 
shock hydrostatic equilibrium sets the early-time light 
curve of our cooling calculations (§11]). 

2.3. Upward propagating detonation 

We compute the evolution during the convective stage 
{th > td) as in WBB. Approximately 30% of the 
is burned to ^^Mg during the convective stage for yb = 
10^2 g cm~2 and a pre-ignition mass fraction X12 = 0.2. 
We then initiate the detonation at the onset of the dy- 
namical runaway [t^ = td) by increasing the tempera- 
ture of the grid point at yb by 1%. In Figure [2] we show 
the evolution of the thermal profile during the convective 
stage as Tb rises {dashed lines) and the hydrodynamic 
stage as the detonation propagates outward [solid lines) . 

The shock wave that defines the head of the detona- 
tion front is sufficiently strong that it triggers complete 
^2C burning in the deepest layers. However, as the shock 
propagates outward into cooler, lower density fuel, it fails 



network solver "pub- 
available online at 



We use F. Timmes' 

lic_aproxl3.f 'j kindly made 

http ://www.cococ ubed.com/code-pages/net_aproxl3.shtml (see 
also ri'imme£|[T99y|T : ' 
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Fig. 2. — Evolution of the thermal profile for yi, = 10^^ g cm~^ 
and an initial composition of 20% ^^C and 80% ^^Fe by mass. The 
sequence of dashed lines show the evolution during the convective 
stage as in WBB. The sequence of solid curves (with times labeled 
in microseconds) show the outward propagating detonation which 
is gradually outrun by the steepening shock. The dotted line shows 
the profile after 9 fis, indicating the maximum extent of the deto- 
nation j/tjet (= 4.3 X 10^^ g cm^^ in this case). 

to trigger ^^C burning on the hydrodynamic timescale 
(~ ^s). The shock ultimately races ahead of the burn- 
ing and the detonation dies at a depth j/dot (see dotted 
line in Figure [2]) . The value of ydct is approximately the 
depth at which the detonation induction length li^^ be- 
comes greater than a scale height h = y/p (this assumes 
the detonation is very nearly planar upon reaching this 
depth; see discussion in S 12. 1|) . 

We now estimate how j/dot depends on the ignition pa- 
rameters. The detonation Mach number M = v/cs « 1.3 
and is close to the Chapman-Jouguet value Mcj = 
(1 + a) + ^1/2^ where a = (7^ - l )Er,uc/2-rhq - 0.1 
7 = 4/3, and i?nuc = 10^'' ergs g"^ ()Landau fc Lifshit j 
(l959t) . Since Und ~ ih, mini the minimum nuclear heating 
time (WBB), we have imd — vt-md ~ Cgth.min, so that 
'ind = h when th.min = ^d- In WBB we derived an es- 
timate of the depth ?/dyn where th,min becomes greater 
than td] thus ydct ~ Udyn, which gives 



J/dot 



(2.4x10" g cm-2) 



0.2 

X12 



3.2 



21 
0.5 



3.0 



2 

514 



(3) 

In Figure [3] we show the composition profile after 
traversal of the shock in the detonation region ydet < 
y < Vb for 2/6,12 = yb/lO^^ g cm"^ = 1. We find that 
for yb,i2 = 0.5 and yh,\-2 = 1, the dominant products 
of burning are the intermediate mass elements ^^Si and 
^^S. For yb,i2 — 5, the dominant products are ^^Si and 
^^S in the top half of the detonation layer's mass and 
^^Ni in the bottom half. As we discuss in § [SJ the in- 
creased crustal heat flux from a crust composed of ^'^Si 
and '^^S can help bring the superburst ignition depth into 
better agreement w ith values inferred from observations 
(jGupta et al.ll2006f ). 



Fig. 3. — Composition profile {solid lines), nuclear energy release 
^-nuc {dashed line), and increase in internal energy A_Eint {dash- 
triple-dot line), after traversal of the shock in the detonation region 
Vdct < y < Vb foi' Vb = lO^'^ g cm~^ and an initial composition 
X12 = 0.2, Xse = 0.8. At depths y < j/dot = 4.3 X lO" g cm-^ the 
detonation fails to trigger ^ C burning and as the shock steepens 
A^int becomes greater than i?nuc. The ^^C at 1; < j/dot will later 
burn in a deflagration (§[4]l. 



We also show in Figure [3] the magnitude of the nu- 
clear energy release in the detonation ii^nuc — 1.1 x 
10^^ ergs g~i and the change in internal energy of the 
material within the detonation region Ai?int — 0.8 x 



10^^ ergs g 1 



Here AE^^ 



int,/ 



E; 



int.O? 



the dif- 



ference between the internal energy, i?i„t./, upon hy- 
drostatic settling after the detonation passes, and the 
internal energy, i^int.Oj at runaway. Thus, a fraction 
1 — (Ai?int/i?nuc) — 0.3 of the nuclear energy release 
is used to power the upward and downward propagat- 
ing shock. The downward propagating shock is weak 
(Ap/p < 0.1) and will penetrate to great depths be- 
fore dissipating. However, the total energy deposited 
< 0.03 MeV nucleon"^ <^ Qcrust, the energy per nucleon 
released in the crust from pycnonuclear and electron cap- 
ture reactions. Thus, the shock is not a significant source 
of crustal heating. 

The boundary at j/dct between the hot ashes and the 
overlying unburned ^^C fuel is convectively unstable. All 
the 12 C between y^ic < y < Vdet will ultimately burn in 
a convective deflagration, perhaps to iron-peak elements 
(§ 2]). Nonetheless, the contrast persists for more than 
a dynamical time and thus does not affect the upward 
propagating shock. 

2.4. Steepening of upward propagating shock 

In Figures|3]and[5]we show the evolution of the temper- 
ature profile and shock overpressure Ap/p = (psh— Po)/Po 
as the shock propagates upward (and downward) . WBB 
showed that as the shock travels down the density gradi- 
ent it steepens as Ap/p cx y~^^^^. Within the detonation 
region (ydot < y < ?/b) we have Ap/p « Ennc/Ei,^t,o ~ 
0.1. This energy release powers an upward and downward 
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Fig. 4. — Evolution of the thermal profile as the shock continues 
to propagate outward (sequence of solid lines, labeled in fis) for 
the same run as Figure [2] The sequence of dashed lines show the 
evolution during the convective stage. The vertical dotted linos 
show the locations of the H/He — ash/C and ash/C — Fe interfaces. 
The dashed-dotted line shows the profile after the shocked layers 
have adiabatically expanded and settled into a new hydrostatic 
solution. The temperature jump due to the shock propagating 
downward into the ^^Fe is too small to be seen on this scale. 

propagating shock such that for y < ydot and y > yb, 



Ap 



2 VS 



int,0 



-9/16 



(X^2 
V 0.2 



1/4 



Vb 



10^^ g cm 2 

9/16 



5/16 



108 g 



(4) 



where the factor of 1/2 in the first expression ac- 
counts for the two shocks. The numerical estimate 
assumes a relativistic degenerate electron gas within 
the detonation region with i?int,o = SYeEp/i = 



(10 ° ergs g" 



-')igim2mYe/0.5), 

117 



where E'f 



is the 

Fermi energy, and i?nuc = (6 x 10^"" ergs g~^)Xi2. Thus, 
Ap/p « 10 upon reaching the H/He layer at t/Ho ~ 
10^ g cm^^. 

When the shock impacts the density discontinuity at 
the H/He — ash/C interface j/hc, a rarefaction wave forms 
and reflects back downward into the previously shocked 
material (see e.g., the curve labeled 8.5 fiS in Figure [5]). 
The upward propagating shock weakens slightly {Ap/p 
decreases by « 20 — 30%) and a steep temperature gradi- 
ent forms at ?/ho- However, despite the steepness of the 
gradient, the fluid is not convectively unstable because 
the interface's initial density contrast is maintained in 
the shock's wake. 

The dashed-dotted line in Figure S] shows the temper- 
ature proflle after passage of the shock and hydrostatic 
readjustment. Since Ap/p is small for y > 10^*^ g cm~^, 
these layers barely need to expand and cool in order to 
regain hydrostatic equilibrium. By contrast, the shock 



Fig. 5. — Evolution of the shock overpressure for the same run as 
the previous figures (we only plot Ap/p > and thus do not show 
the rarefaction wave that trails the shock). The shock steepens as 
roughly y~^f^^ . The vertical dotted lines show the locations of the 
H/He — ash/C and ash/C — Fe interfaces. 

is strong for y < 10^ g cm~^ and these layers expand 
and cool significantly. In § 2] we show that the cooling 
of the region 10^ g cm~^ ^ y ^ 10^ g cm~^ results in 
a brief (w 0.1 s), bright, flash that precedes the precur- 
sor burst by several seconds. The duration of the flash 
is determined by the hydrostatic temperature profile in 
this region, which we find is nearly isothermal for the 
following reason. The layer is initially radiative Tq cx 



and the shock is gas pressure-dominated. 



1/4 1/4 

Tsh/Ta (xpsh/po- Thus, equation (g]) gives Tsh oc y"^/^'^. 
The layer then adiabatically expands from psh back to 
P = Pa ^ gy (the layer always remains plane-parallel) and 
the temperature decreases to a hydrostatic value that 
varies weakly with depth Thso ~ Tsh{po/Psh)'^^^ y~'^^^^. 

Near the top of our grid, the shock becomes radia- 
tion pressure dominated with a post-shock temperature 
Tsh ~ (SAp/ay^^ (X y'^/S''. The diffusion time of the 
photons generated in the shock t^ie nyh/c is larger 
than the shock propagation time igh h/vsh as long as 
y ^ c/nvsh, where the opacity k ~ 0.2 cm^ g~^. Thus, 
photons do not preheat the gas ahead of the shock until 
y ^ 10^ g cm^^ and radiative transport can be neglected 
over the range of our shock calculation. 

For a strong shock, hP = (fsh/cs)^ ~ Ap/p. In 

the radiative outer layers Cg « (10^ cm s~^)yg^* so that 

Wsh ~ (10^ cm s"^)% where y^ = y/lO^ g cm"^ 
and we used ([4]). Since the neutron star escape speed 
Ucsc ~ c/3, the shock ejects only the outermost layers 
?/ < 1 g cm~^. 

3. HYDROGEN AND HELIUM BURNING 

As the shock propagates through the upper layers, it 
deposits considerable entropy. In § 13.11 we show that if 
the superburst occurs during the latter half of the Type I 
bursting cycle, there is enough ^He present in the H/He 
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layer that the shock triggers unstable ^He burning. The 
^He is consumed within only a few seconds and, as we 
show in § 14.21 this rapid energy release results in a nearly 
Eddington-limited Type I burst. In ^ 13.21 we show that 
if the superburst occurs during the first half of the cycle, 
the ^He still burns, but only in a thermally stable man- 
ner, over a duration of > 100 s. Such stable ^He burning 
does not result in a Type I burst because the ^He is con- 
sumed so slowly that the heat flux from burning never 
gets very high (though we show that stable rp-process H 
burning may increase the flux considerably). 

3.1. Shock- Triggered Unstable Helium Burning 

The shock will trigger unstable '^He burning only if 
the base of the ^He layer is thermally unstable after 
settling into its post-shock hydrostatic configuration. 
This is because the shocked layers hydrostatically settle 
within only a few dynamical times 10/is). By con- 
trast, the heating timescale from unstable ^He burning 
at these depths is al ways greater than ^ 1 ms (see e.g., 
[Weinberg et~ani2006bl . 

In Figure [6] we show the hydrostatic post-shock ^He 
ignition curve as a function of Ap/p at j/He- We de- 
fine ignition according to the one-zone thermal insta- 
bility criterio n de^g/dT = dccooi/ dT ([ Fuiimoto ct al.l 
[1981: Fus hiki fc LambllOSTtlCumniing fc B ildstcn 200Q) , 
where e^a is the triple-alpha energy generation rate (we 
use that from lFushiki fc Lamblll987f ) and ecooi = pKT/y"^ 
is an approximation to the local cooling rate for a ther- 
mal conductivity K. To the right of the ignition curve 
desa/dT > dtcoo\/dT and burning is unstable. The 
three curves are for different local accretion rates per 
unit area rh, in units of the local Eddington rate rriEdd- 
The lower the m, the larger the ^He abundance at t/hc 
(jCumming fc BildstenI l200df ) , and hence the lower the 
Ap/p threshold for ignition. If the last Type I burst was 
too recent (i.e., the phase, and thus j/He, is too small), 
there is no value of Ap/p capable of triggering unsta- 
ble burning si nce at low densities c ooling is always faster 
than heating (jFuiimoto et aLlfTOBlh . If Ap/p is too large, 
^He burning is always stable and will not produce a Type 
I burst; however, such a strong shock is unlikely since it 
requires j/f, > 5 x 10^^ g cm~^ and a superburst early in 
the Type I burst cycle. 

For a ^^C detonation at j/b > 5x10^"'^ g cm~^, the shock 
is sufficiently strong to trigger unstable ^He burning for 
the latter « 40% — 50% of the burst cycle. If a less 
massive layer of "^He is present, a Type I burst will not 
be triggered. 

A precursor has been seen in each of the five cases 
where the onset of the burst was detected, and it is un- 
likely that all are due to shock-triggered '*He burning. 
However, as we describe in § HI some of the observed 
precursors may not be due to unstable ^He burning but 
rather to ^^C burning at depths just below ?/hc- Further- 
more, we may be underestimating the fraction of super- 
bursts that trigger unstable "^He burning for two reasons. 
First, the H/He layer may be more He-rich than our cal- 
culations assume, especially at higher m. This is be- 
cause superbursts tend to occur in systems with large 
values of a (> 1000), the rat io of the time-ayerage d 
persistent to burst luminosity ijin't Zand et al.l l2004af ) . 
The Type I bursts in these systems also tend to be of 
short-duration, consistent with a pure "^He flash. On the 
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Fig. 6.— The "He ignition curve (solid lines) as a function of 
the shock overpressure Ap/p at yjie for '"/"lEdd = 0.05,0.1, and 
0.3. The dashed hnes indicate the overpressure at ^ho for a- ^'^C 
detonation ignited at yb,i2 = 0.5, 1, and 5 with X12 = 0.2. The 
two squares along each curve indicate a Type I burst cycle phase 
of = 0.6 and (p = 0.8. The squares on the 0.05mEdd curve 
are labeled by the temperature T/10*K at j;hc after hydrostatic 
settling. 

other hand, the large a values may also indicate that 
most of the ^He is burning stably ( Naravan & Hcyl 200^; 
ICooper fc Naravanll2005l : ICooper et all 120061. and it is 
not clear how a shock would trigger a transition from sta- 
ble to unstable burning. Second, if a convective ^^C de- 
flagration propagates outward from j/dct and reaches yuc , 
it may also trigger unstable *Hc burning. Such a defla- 
gration would propagate at roughly the turbu lent speed 
vturh > 10^ cm s-1 (cf. iTimmes fc NiemeveijfeOOO ) and 
reach the "^He layer a few milliseconds after the shock. 

3.2. Stable Hydrogen and Helium Burning 

If the shock does not trigger unstable ''He burning, the 
H and ^He still burn in a stable manner due to the large 
heat flux from the ashes of ^^C burning below. We ac- 
count for three sources of stable burning: hot CNO burn- 
ing chcno — 5.8 X lO'^ZcNO ergs g"-' (we assume 
the initial mass fraction of CNO nuclei is .Z'cno — 0.01), 
triple-alpha '*He burning 63^, and H burning via the rapid 
proton capture (rp) process erp of [Wallace fc WooslevI 
fl98l). The evolution of stable H and ^He burning dur- 
ing the flrst 1000 s following hydrostatic settling is shown 
in the top panel of Figure [T] 

To calculate Cip, note that after hydrostatic settling 
the temperature and density at the base of the H/He 
layer ?/hc are T ~ 5 x 10^ K and p ~ 10*^ g cm"^. The 
primary breakout reactio n from the hot CN O cycle is 
therefore ^^0{a,-fy^Ne (jSchatz et al.lfl999h .^ As the 

* Due to the heat flux from the ashes of ^^C burning, the tem- 
perature at j/Hc reaches T ~ 9 X 10* K just 0.01 s after hydrostatic 
settling. At this t emperature, a secon d hot CNO cycle emerges 
that includes ^*Ne llWiescher et al.lll999l) . and the primary break- 
out reaction may instead be ^*Ne(a, pl'^^Na JSchatz et al.|[T999l) . 
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Fig. 7. — Stable H and ^He burning at the base of the H/He layer 
j/He as a function of time after hydrostatic readjustment for j/t, 12 = 
1, m = O.lriiEddi S'lid a- Type I burst phase = 0.5 (i.e., the shock 
does not trigger unstable "^He burning). The top panel shows the 
energy generation rate {solid lines; in units of 10"*^** ergs g"-*^ 
due to rp-process burning Crp (equation [5]), stable triple-alpha *He 
burning e^a, and hot CNO burning ehcno- The dashed lines show 
the mass fractions Xjj, ^Hci and ^CNO- Shown for comparison 
(dash-dot lines) are ega and Xjje during the convective stage of 
unstable ''He burning for the case yi, 12 = li = O.lrriEddi and 
(f> = 1. The bottom panel shows the temperature Tg = K 
at 3/Hei and the luminosity L^g = ergs s~^ at the surface. 

The set of three Tg and L38 curves show, from bottom to top, the 
result if in equation ^ we take Enuc = 0, e^uc = esa + eHCNOi and 

Enuc = Erp + eSa + tHCNO- 

thermal wave from the hot ashes of ^^C burning dif- 
fuses through the H/He layer, the temperature at yuc 
rapidly rises and the -'^^Ne captures a proton before it 
can P decay {ti/2 = 17.2 s) and return to the hot CNO 
cycle. A series of fast proton captures ensue until the nu- 
clear flow reaches the first waiting point of the rp-process, 
24si = 0.102 s: IWiescher et all 11998). Energy pro- 
duction by rp-processing drops until higher temperatures 
initiate further proton and a captures. Here we follow 



Wiescher et al. (1999; see also lCooper fc Naravan|[2006l) 
and approximate the energy generation rate of rp-process 
burning by accounting for just the energy release in burn- 
ing to 2*Si. Since ^^0(a,7)^^Ne is the slowest reaction 
of the sequence, it governs the total reaction rate of the 
flow, and 

e,p48 = (0.018)XHe^cNoQrppr9"'^'e-5«5/^«, (5) 

in units of 10^^ ergs g^^ s"^ Here Qrp = 29.96 MeV 
is the energy release in burning to ^"^Si and the 
^^0(a,7)^^Ne reaction rate is from Caughla n fc Fowleil 
([l98_8,). This approximation only allows us to place a 
lower-bound on the contribution of rp-process burning to 
the light curve; the full calculation is left for future study. 
Nonetheless, as the bottom panel of Figure [7] shows, it 
has a noticeable impact on the early time light curve 
and if the shock does not trigger unstable ''He burning, 
Crp/esa ~ 5 — 30 over the duration of the burst. 

4. COOLING MODELS AND LIGHT CURVES 

Following the passage of the shock and hydrostatic 
readjustment, a cooling wave propagates inward from the 
surface and the shock-heated layers begin to cool via heat 
diffusion, 

dT _dF 

where the flux F = pK{dT/dy), Cnuc is the nuclear en- 
ergy generation rate due to any remaining fuel, and ei, is 
the neutrino energy loss rate. The light curve is set by 
the final burn temperature profile Tj (y) , input as the ini- 
tial condition of ([6]), and the heat fiux from any burning 
Cnuc following hydrostatic settling. 
Previous studies assumed Tf{y) was given by 

the total energy release Jj,^ CpdT — Xi2E^^c — 

(lO'^^ergs g^^)Ari2 assuming an instantaneous and com- 
plete isoba ric ^^C b urn between 10^ g c m"^ < y < yj,. 
This gives (jCumming fc Macbeth, 2004; .Gumming et all 
[2001), 

T,,.oba.-(5.1xlO^K)(|f)'''(^)'''. (7) 

Here we instead use the results of our time-dependent 
detonation calculation to determine Tf{y). 

The Cnuc term has contributions from two regions of 
unburned fuel: the ^^C between ^hc < y < 2/det and the 
H/He at y < i/hc- As noted in § 12.31 due to the sharp 
temperature gradient at j/dot, the ^^C will likely burn 
in a convective defiagration that will reach the H/He 
layer a few ms after passage of the shock. For simplicity, 
we assume an instantaneous isobaric ^^C burn (approx- 
imately given by equation [7] ) in the deflagration region 
J/Hc < y < ydet- We assume an energy release in the 
deflagration of Enuc = 10^'' ergs g~^, corresponding to 
^^C burning to ^^Fe with X12 ~ 0.1, the fraction of fuel 
remaining at the onset of the runaway that follows the 
convective stage. 

To account for H and "^He burning, we adopt separate 
strategies depending on whether the shock triggers un- 
stable "^He burning. If the shock triggers unstable '^He 
burning f§ 13. ip . the burning layer quickly becomes con- 
vective and its thermal evolution during the rise is not 
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Fig. 8. — Temperature profiles for y^, = 10 g cm , rh = O.lriiEddi a^nd a Type I burst phase = 0.5. The mass of freshly accreted 
H/He is too small for the shock to trigger unstable *He burning (see Figure |6ll. Shown are profiles at ignition (dashed- dotted line), after 
the '^■^C detonation and shock breakout but before the isobaric deflagration (dashed line), and at ten times after hydrostatic settling: 
t = 0,0.005,0.05,0.3, 1, 10, 1000, 10*, 6 x 10^,2 x 10^ s. The left panel shows the cooling profiles over the entire range in y and the right 
panel is a zoom in of the H/He layer with the ten profiles labeled in sequential order. The vertical dotted lines show the locations of the 
H/He — ash and ash — Fe interfaces. 



described by the heat diffusion equation Instead, 
we solve the rise as in Weinberg et al. (2006b; here we 
account only for a-capture reactions) , which we then su- 
perimpose onto the cooling solution. If the shock does 
not trigger unstable ^He burning, the H and "^He still 
burn in a stable manner (§ l3.2p . In this case a convective 
zone does not form and the thermal evolution is described 
by equation We set Cnuc = ehcno + esa + £rp , the 
three sources of stable burning described in § 13.21 and 
solve for the rate of change of the composition dXi (y) /dt 
simultaneously with the diffusion equation (see Figure 
III). 

We use the method of lines to solve equation 
setting dlnT/dlny = 1/4 at the outer boundary 
(10** g cm~^) and dT/dy — (vanishing flux) at the inner 
boundary (10^'' g cm^^). In Figures [8] and [9] we show the 
cooling for yb,i2 = 1, rh = O.lmEdd, and a Type I burst 
phase (/) = 0.5 and 0=1, respectively. In the former 
case, the shock does not trigger unstable ^He burning. 
The corresponding light curves are shown in Figures 1101 
and[TlJ We find that in all cases the light curve evolves 
through three stages: a shock breakout stage {t < 1 s; 
§ 14. ip . a precursor stage (1 s < < < 20 s; § I4.2p . and a 
superburst stage {t > 20 s; S I4.3p . 

4.1. Shock breakout 

For times < 1 s after the detonation, the cooling of the 
outermost shock-heated layers {y < 10^ g cm~^) domi- 
nates the emission. It produces a brief (« 0.1 s), bright, 
flash that precedes the precursor burst peak by « 3 — 10 s 
and has an initial luminosity that is super-Eddington. 
The evolution of the temperature profile within the H/He 
layer during this cooling epoch is shown in the right panel 
of Figure [8] (lines {0, 1, 2}). The luminosity of the fiash 
decays as a power- law with time L oc i^^/^, which can 



be understood as follows (| Gumming et al.l |2006| give a 
similar explanation for the superburst cooling at times 
100 s < t < 1000 s). The cooling wave propagates in- 
ward from the surface and reaches a column depth at 
time t w ithcrm(j/*), the thermal timescale of the layer. 
For y < y* the atmosphere has an almost constant heat 
flux with depth and an approximately constant opac- 
ity K so that L oc T^/j/ = Tj{y,,)/y^, — constant. 
For y > y^, the temperature profile is nearly unchanged 
from its initial state Tf{y). Since initially the region 
10^ g cm~^ < y < J/He is nearly isothermal (see § 12. 4p . 
L oc y^^, and since gas pressure dominates tthcrm — 
3CpKy^/4acT^ (xy^,so that L oc t''^^'^. 

As the Type I burst phase (p ^ 1, the delay between 
the onset of the flash and the peak of the precursor burst 
increases slightly (Figure [TOl) because the larger j/hc is, 
the greater ithorm(yHc) is. As the ^^G ignition depth j/f, 
increases, the fiash becomes more luminous (Figure llip 
because the larger yh is the larger Ap/p is at a given 
depth y (equation [4]), and thus the hotter the envelope 
is after hydrostatic settling. 

While the outermost layers y < 10^ g cm~^ cool, the 
heat fiux from the hot ashes of ^^G burning begins to 
diffuse outward and heat up the layer 10^ g cm~^ ^ 2/ < 
j/Hc- A steep temperature gradient develops near j/He 
and a convective zone forms. Since we assume the energy 
transport is described by heat diffusion (equation [6]), we 
underestimate the efficiency of energy transport during 
this brief (< 1 s) convective epoch. The delay between 
the peak of the shock breakout fiash and the peak of the 
precursor burst may therefore be somewhat shorter than 
we have found and the valley between the peaks may not 
be as deep. 

Unlike the unstable ^He burning precursor, which 
might be triggered by a shock or an incident ^^G defia- 
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Fig. 9. — Temperature profiles for yi, = lO^'^ g cm~^, m = 
O.lrhEddi and a Type I burst phase = 1. Here the shock 
triggers unstable ^He burning, unlike in Figure |8] Shown are 
profiles at several times after hydrostatic settling: Lines {0, 1, 2} 
show the prompt cooling just after shock breakout at times 
{0,0.01,0.1 s}. Lines {3,4,5,6,7} show the H/He layer heating 
up during ^He burning: label 3 shows the growing convective re- 
gion at {0.6, 0.9, 1.2, 1.4, 1.7 s}, and Unes {4, 5, 6, 7} show the burn- 
ing as the radiative layers heat up and the fiux begins to rise at 
{3.1,3.9,4.5,5.3 s}. Lines {8,9,10} (dash-triple-dot lines) show 
the superburst cooling after "^He exhaustion at {15, 1000, 3x 10** s}. 
The dashed-dotted line is the profile at ignition and the dashed line 
the post-shock, pre-deflagration, profile. 

gration, the flash is perhaps a unique signature of a "'^^C 
detonation. A sub-second spike is seen at the onset of 
the superburst in both cas es in which RX TE caught the 
rise (4U 1820-30: .Strohmaver fc Bro wj [2002,; 4U 1636- 
54: iStrohmaver fc Markwardtll2002f) . In both cases, the 
peak of the spike precedes the peak of the precursor by 
~ 2 — 3 s. Although consistent with the features of a 
shock breakout, it is possible the spike is instead just 
the first peak of a photospheric radius expansion burst. 
However, in the case of 4U 1636-54, the peak flux of 
the precursor is approximately 60% fainter th an that of 
the b rightest Type I bursts from this system (jKuulkersI 
120041) . This suggests that at least in 4U 1636-54, the 
precursor is not Eddington-limited and thus not a radius 
expansion burst. 

4.2. Precursor hurst 

For 1 s < t < 20 s, the rise and decay of the unsta- 
ble ^He burning precursor, and the cooling of ashes of 
^^C burning just below t/ho, dominate the emission. If 
unstable '*He burning is triggered, it results in a Type 
I burst that begins to rise « 1 s after the peak of the 
shock breakout flash. During the first « 1 s of ^He burn- 
ing, a convective zone develops and moves outward from 
j/Hc to lower pressure (line 3 in Figure [9]). The energy 
generation rate £30, is so high during the initial stages 
of burning that the timescale for the convective zone to 
grow is much shorter than the thermal timescale at the 
convective-radiative interface. As a result, the outer ra- 
diative layers do not get heated during the first second 
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Fig. 10. — Luminosity as a function of time during the early 
stages of the cooling for an ignition depth Vb 12 = rn = O.lrhEddi 
and a Type I burst phase <p = 0.5 (dash-triple- dot), rj> = 0.8 
(dashed), and = 1 (solid). For = 0.8 and = 1 the shock 
triggers unstable *He burning. For = 0.5, H and ^He are burned 
only in a stable manner over 100 s; approximately 90% of the 
heat fiux at peak is due to the cooling ashes of ^^C burning at 
depths just below yne (see Figure [Tjl. 

of ^He burning and the cooling wave continues to prop- 
agate inwards through the outer shock-heated layers. A 
steep temperature gradient develops at the convective- 
radiative interface due to the large compositional con- 
trast between the ^He-ric h matter that i s burning and 
the outer H-rich material (iWeinberg et al.l , 2006b) . Even- 
tually, the thermal time at the interface becomes shorter 
than the convective growth time and the convective zone 
slowly retreats back to the base yne (lines {4, 5, 6, 7} in 
Figure [9]). The radiative layers then finally start to heat 
up and the flux begins to rise, reaching a peak luminosity 
after several seconds. 

The peak luminosity of the Type I burst increases as 
0^1 and the amount of ^He available to burn increases 
(Figure fTO|) . After the peak, the H/He layers cool and 
the thermal wave continues its inward march, gradually 
penetrating the ashes of ^^C burning (lines {8,9,10} in 
Figure [91). 

Importantly, even if unstable ^He burning is not trig- 
gered, the cooling ashes of ^^C burning just below j/hc 
results in a flux rise and decay that is similar in appear- 
ance to that of a Type I burst, though with a peak flux 
that is smaller by a factor of w 2. This can be seen 
in the = 0.5 curve of Figure [TOl a case in which the 
shock fails to trigger unstable "^He burning (also com- 
pare lines {3, 4, 5} of Figure [5] with lines {3, 4, 5, 6, 7} of 
Figure [9]). In this case, the Type I- like rise and decay 
is powered by the hot ashes of ^^C burning, although 
stable nuclear burning of H and ^He also contribute, in- 
creasing the peak flux by w 10% — 20% (see the bottom 
panel of Figure [Tj this may be an underestimate since 
we only account for stable H burning via the rp-process 
up to ^*Si). Thus, even if there is no '^He layer, if the 

C at y « 10^ g cmT^ burns there will still be a pre- 



10 



WEINBERG & BILDSTEN 



O 



O 

.S 0.1 



3 



shock breakout + precursor 



superburst 



10 




1 10 100 

Time (seconds) 



1000 



Fig. 11. — Luminosity as a function of time for ignition depths j/t,,i2 = 0.5 (dash-triple-dot), 1 (solid), and 5 (dashed), rh = O.lrh^^^, 
and a Type I burst phase = 1. In each case, the shock triggers unstable ^He burning and the peak of the shock breakout precedes the 
peak of the precursor burst by a few seconds. The vertical dotted line demarcates where the light curve transitions from the shock breakout 
and precursor/Type I burst to the superburst. 



cursor burst, albeit one that is about half as bright as an 
ordinary Type I burst. 

In three of the five cases in which a precursor was 
seen, the peak flux of the precursor was smaller than 
the brightest of the system's ordinary bursts by a fac- 
tor of 1.5 — 2 (the exceptions are the precursor from 
4U 1820-30, which looked like an ordinary burst, and 
4U 1254-69 which was brighter than an ordinary burst; 
lKuulkerg|[200llCumming et al.ll2006f ). It is thus possible 
that some of the observed precursors are not the result 
of unstable "^He burning. 

4.3. Late-time light curve 

For t > 20 s, the light curve is powered by the con- 
tinued inward propagation of the cooling wave into the 
ashes beneath yue (se e Figu re [H lines {6,7,8,9}). As 
iCumming fc MacbethI (|2004[ ) showed, the light curve de- 
cays as a broken power-law, with the break marking the 
time when the cooling wave first reaches the deeper 
yb is, the later the break (Figure [TT|) . Our final tempera- 
ture profile Tf{y) in the detonation region j/det < y < yt 
is only slightly steeper than the profile in the isobaric 
deflagration region ?/He < D < Udet- Thus, our computed 
light curves at times t > 10 s are similar to those of 
ICumming fc MacbethI (|2004f ). who assumed an isobaric 
burn throughout. 

Here, as in previous studies, the calculated slope of the 
pre-break light curve {t « 100 — 1000 s) is steeper than 
the o bserved slope (see e.g., Figure 6 of ICumming et all 
l2006l) . The emission during this epoch is powered by the 
cooling of the deflagration- heated region j/ho < y < J/det, 
and in order to better match the observations, Tf{y) in 
this region must be steeper than we have assumed. In- 
deed, since the deflagration is likely to be highly turbu- 
lent due to the strong gravity and the hour-long convec- 
tive churning that precedes the runaway, our assumption 
of a complete isobaric burn may not be correct. The 



deflagration may instead leave behind increasingly large 
pockets of unburned fuel as it accelerates down the den- 
sity gradient, thereby yielding the steepened temperature 
gradient implied by the observed light curve. 

5. SUMMARY AND CONCLUSIONS 

We modeled the hydrodynamic ^^C combustion wave 
that forms during a superburst as a detonation. We 
found that the detonation propagates through the deep- 
est layers and drives a shock wave that steepens as it 
travels upward. Upon reaching the H/He layer, the shock 
is sufficiently strong that it triggers unstable ^He burning 
if the superburst occurs during the latter half of the Type 
I burst cycle. We showed that the light curve that results 
from a shock-triggered "^He burn is similar to that of a 
normal Type I burst. The main difference is that the 
peak luminosity is somewhat smaller than an ordinary 
burst's if the superburst occurs early in the second half 
of the Type I cycle, when there is less ''He present. 

A precursor burst has been found in each of the five 
cases where the onset of the superburst was observed. 
Given our results, it is unlikely that all are due to shock- 
triggered unstable ^He burning. It is possible that our 
analysis underestimates the fraction of the cycle over 
which the shock ignites the ''He (see discussion in § [3]) . 
However, because the short thermal time at low densities 
results in a sharp turnover in the ^He ignition curve at 
j/Hc < 10^ g cm~^, unstable ''He ignition is unlikely over 
much more than half of the cycle, even in a layer of pure 
^He. A more likely explanation is that some fraction of 
the observed precursor bursts are due to ^^C burning at 
depths just below the H/He layer (rp-process H burning 
may also be stronger than we have found). We showed 
that the light curve of a burst in which unstable ''He 
burning is not triggered still has a Type I-like precursor, 
but with a peak luminosity that is smaller by a factor of 
w 2. This may explain why the peak luminosity of three 
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of the five observed precursors was about half as bright 
as the system's brightest ordinary bursts. A larger sam- 
ple of precursor bursts is needed in order to better test 
this hypothesis. 

Although we modeled the combustion wave as a deto- 
nation, the entire ^^C layer may instead burn in a defla- 
gration (see ? 12. ip . How might observations distinguish 
between the two modes of propagation? The late-time 
light curve, though well observed, is unlikely to help since 
at late times the cooling of the ashes of a detonation ap- 
pears very similar to that of a deflagration f§ I4.3p . The 
bright precursor bursts also do not imply a detonation, as 
a ^^C deflagration may also trigger unstable ''He burning 
when it impacts the H/Hc layer. Thus, perhaps the most 
promising signature of a detonation is the shock break- 
out, a bright, « 0.1 s, flash that precedes the precursor 
burst by ~ 3 — 10 s. This may be the origin of the spike 
seen at the onset of the burst rise in 4U 1820-30 and 4U 
1636-54. 

The superburst ignition depth is sen sitive to the 
heat release from deep crustal reactions (jBrownl 120041 : 
ICumming et al.l l2006l ). The latter, in turn, de- 
pends on the ash composition of previous superbursts 
(|Gupta et al.ll200^ . Using a 13 isotope a-chain network, 
we found that the ashes of the detonation consist primar- 
ily of ^^Si, •^^S, and '^^Ar . This is potentially interesting 
since iGupta et al.l (|2006f ) find that a crust composed of 
nuclei with A = 32 — 44, as compared to A ~ 50 — 60, 
has a larger heat flux from electron captures. It results 



in a significantly hotter crust and decreases the super- 
burst ignition depth yh by a factor of « 2, which helps 
bring yb into better agreement with values inferred from 
observations. Further study is required as our network 
is small and the results likely depend on the initial com- 
position of the ignition layer (we assumed a ■^'^C/^^Fe 
mixture). For example, superbursts may ignite in a sea 
of light nuclei if most of the ^^C that fuels the super- 
burst is produced in stable ''He b urning (jNaravan fc Hevl 
I2003t ICooper fc Naraviml |2005[ ) . Alternatively, super- 
bursts may ignite in a sea of very hea vy nuclei (A ^ 100) 
made during rp-process H burning (ISchatz et al.l [19991 . 
I2001h . In the latter case. ISchatz et all (| 2003f) found that 
the superburst nuclear energy release is dominated by the 
photodisintegration of the heavy nuclei, which results in 
a final ash composition that instead consists primarily of 
iron-group elements. 
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